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ABSTRACT 

We present spectra for 33 previously unclassified white dwarf systems brighter 
than V = 17 primarily in the southern hemisphere. Of these new systems, 26 are 
DA, 4 are DC, 2 are DZ, and 1 is DQ. We suspect three of these systems are unre- 
solved double degenerates. We obtained VRI photometry for these 33 objects as 
well as for 23 known white dwarf systems without trigonometric parallaxes, also 
primarily in the southern hemisphere. For the 56 objects, we converted the pho- 
tometry values to fiuxes and fit them to a spectral energy distribution using the 
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spectroscopy to determine which model to use (i.e. pure hydrogen, pure hehum, 
or metal-rich helium), resulting in estimates of T^g and distance. Eight of the 
new and 12 known systems are estimated to be within the NStars and Catalogue 
of Nearby Stars (CNS) horizons of 25 pc, constituting a potential 18% increase 
in the nearby white dwarf sample. Trigonometric parallax determinations are 
underway via CTIOPI for these 20 systems. 

One of the DCs is cool so that it displays absorption in the near infrared. 
Using the distance determined via trigonometric parallax, we are able to constrain 
the model-dependent physical parameters and find that this object is most likely a 
mixed H/He atmosphere white dwarf similar to other cool white dwarfs identified 
in recent years with significant absorption in the infrared due to collision-induced 
absorptions by molecular hydrogen. 

Subject headings: solar neighborhood — white dwarfs — stars: evolution — 
stars: distances — stars: statistics 



1. Introduction 

The study of white dwarfs (WDs) provides insight to understanding WD formation 
rates, evolution, and space density. Cool WDs, in particular, provide limits on the age of the 
Galactic disk and could represent some unknown fraction of the Galactic halo dark matter. 
Individually, nearby WDs are excellent candidates for astrometric planetary searches because 
the astrometric signature is greater than for an identical WD system more distant. As a 
population, a complete volume limited sample is necessary to provide unbiased statistics; 
however, their intrinsic faintness has allowed some to escape detection. 

Of the 18 WDs with trigonometric parallaxes placing them within 10 pc of the Sun (the 
RECONS sample), all but one have proper motions greater than I'.'O yr~^ (94%). By com- 
parison, of the 230 main sequence systems (as of 01 January 2007) in the RECONS sample, 
50% have proper motions greater than I'/O yr~^. We have begun an effort to reduce this 
apparent selection bias against slower-moving WDs to complete the census of nearby WDs. 
This effort includes spectroscopic, photometric, and astrometric initiatives to characterize 
newly discovered as well as known WDs without trigonometric parallaxes. Utilizing the Su- 
perCOSMOS Sky Survey (SSS) for plate magnitude and proper motion information coupled 
with data from other recently published proper motion surveys (primarily in the southern 
hemisphere), we have identified relatively bright WD candidates via reduced proper motion 
diagrams. 
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In this paper, we present spectra for 33 newly discovered WD systems brighter than 
V = 17.0. Once an object is spectroscopically confirmed to be a WD (in this paper for 
the first time or elsewhere in the literature), we obtain CCD photometry to derive Tgfj and 
estimate its distance using a spectral energy distribution (SED ) fit and a model a tmosphere 
analysis. If an object's distance estimate is within the NStars (IHenry et al.ll2003l ) and CNS 
( iGliese fc Jahreifilll99ll ) horizons of 25 pc, it is then added to CTIOPI (Cerro To lolo Inter 



American Observatory P arallax Investigation) to determine its true distance (e.g. |Jao et al. 



20051 . iHenrv et al.ll2006f ). 



2. Candidate Selection 



W e used recent high pr oper motion (HPM) surveys ( Pokorny et al.ll2004j : ISubasavage et al 
2005al lbl: I Finch et al.l 120071 ) in the southern hemisphere for this work because our long-term 
astrometric observing program CTIOPI, is based in Chile. To select good WD candidates 
for spectroscopic observations, plate magnitudes via SSS and 2MASS JHKs are extracted 
for HPM objects. Each object's (-R59F — J) color and reduced proper motion (RPM) are 
then plotted. RPM correlates proper motion with proximity, which is certainly not always 
true; however, it is effective at separating WDs from subdwarfs and main sequence stars. 
Figure [T] displays an RPM diagram for the 33 new WDs presented here. To serve as examples 
for the locations of subdwarfs and main sequence stars, recent HPM dis coveries from the 



Super COSMOS-RECONS (SCR) proper motion survey are also plotted (ISubasavage et al. 



2005al lb[). The solid line represents a somewhat arbitrary cutoff separating subdwarfs and 
WDs. Targets are selected from the region below the solid line. Note there are four stars 
below this line that are not represented with asterisks. Three have recently been spectro- 
scopically confirmed as WDs (Subasavage et al., in preparation) and one as a subdwarf (SCR 
1227- -4541, denoted by "sd") that fell just below the line at (i?59F - J) = 1-4 and Hr^^^ = 
19.8 JSubasavage et al.ll2005b[ ). 



Completeness limits (S/N > 1 0) for 2MASS are J = 15.8, H = 15.1, and Ks = 14.3 
for uncontaminated point sources (jSkrutskie et al.ll2006l ). The use of J provides a more 
reliable RPM diagram color for objects more than a magnitude fainter than the Ks limit, 
which is particularly important for the WDs (with (J — Ks) < 0.4) discussed here. Only 
objects bright enough to have 2MASS magnitudes are included in Figure [H Consequently, 
all WD candidates are brighter than V ~ 17, and are therefore likely to be nearby. Objects 
that fall in the WD region of the RPM diagram were cross-referenced with SIMBAD and 
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McCook fc Sio"3 (|l93)0 to determine those that were previously classified as WDs. The 
remainder were targeted for spectroscopic confirmation. 

The remaining 33 candidates comprise the "new sample" whose spectra are presented 
in this work, while the "known sample" constitutes the 23 previously identified WD systems 
without trigonometric parallaxes for which we have complete VRIJHKs data. 



3. Data and Observations 

3.1. Astrometry and Nomenclature 

The traditional naming convention for WDs uses the object's epoch 1950 equinox 1950 
coordinates. Coordinates for the new sample were extracted from 2MASS along with the 
Julian date of observation. These coordinates were adjusted to account for proper motion 
from the epoch of 2MASS observation to epoch 2000 (hence epoch 2000 equinox 2000). 
The coordinates were then transformed to equinox 1950 coordinates using the IRAF proce- 
dure precess. Finally, the coordinates were again adjusted (opposite the direction of proper 
motion) to obtain epoch 1950 equinox 1950 coordinates. 

Proper motions were taken from various proper motion surveys in addition to unpub- 
lished values obtained via the SCR proper motion survey while recovering previously known 
HPM objects. Appendix A contains the proper motions used for coordinate sliding as well 
as J2000 coordinates and alternate names. 



3.2. Spectroscopy 

Spectroscopic observations were taken on five separate observing runs in 2003 Octo- 
ber and December, 2004 March and September, and 2006 May at the Cerro Tololo Inter- 
American Observatory (CTIO) 1.5m telescope as part of the SMARTS Consortium. The 
Ritchey-Chretien Spectrograph and Loral 1200x800 CCD detector were used with grating 
09, providing 8.6 A resolution and wavelength coverage from 3500 to 6900 A. Observations 
consisted of two exposures (typically 20 - 30 minutes each) to permit cosmic ray rejection, 
followed by a comparison HeAr lamp exposure to calibrate wavelength for each object. Bias 
subtraction, dome/sky fiat-fielding, and extraction of spectra were performed using standard 
IRAF packages. 



The current web based catalog can be found at 



http://heasarc.nasa.gov/W3Browse/all/mcksion.html 



- 5 - 



A slit width of 2" was used for the 2003 and 2004 observing runs. Some of these data have 
flux calibration problems because the slit was not rotated to be aligned along the direction 
of atmospheric refraction. In conjunction with telescope "jitter", light was sometimes lost 
preferentially at the red end or the blue end for these data. 

A slit width of 6", used for the 2006 May run, eliminated most of the flux calibration 
problems even though the slit was not rotated. All observations were taken at an airmass 
of less than 2.0. Withi n our wavelength window, the maximum atmospheric differential 



refraction is less than 3" ( lFilippenkdll982l ) . A test was performed to verify that no resolution 
was lost by taking spectra of a F dwarf with sharp absorption lines from slit widths of 2" to 
10" in 2" increments. Indeed, no resolution was lost. 

Spectra for the new DA WDs with T^s > 10000 K are plotted in Figure [2] while spectra 
for the new DA WDs with T^s < 10000 K are plotted in Figure [31 Featureless DC spectra 
are plotted in Figure HI Spectral plots as well as model fits for unusual objects are described 
in 5 4.2. 



3.3. Photometry 



Optical VRI (Johnson V, Kron-Cousins RI) for the new and known samples was ob- 
tained using the CTIO 0.9 m telescope during several observing runs from 2003 through 
2006 as part of the Small and Moderate Aperture Research Telescope System (SMARTS) 
Consortium. The 2048x 2046 Tektronix CCD caniera w as u sed with the Te k 2 VRI filter 
setl^. Standard stars from Graham ( 1982 ). Bessel (1990), and Landolt ( 1992 ) were observed 
each night through a range of airmasses to calibrate fluxes to the Johnson-Kron-Cousins 
system and to calculate extinction corrections. 

Bias subtraction and dome fiat-fielding (using calibration frames taken at the beginning 
of each night) were performed using standard IRAF packages. When possible, an aperture 
14" in d iameter was us ed to determine the stellar flux, which is consistent with the aperture 
used by iLandoltl (119921 ) for the standard stars. If cosmic rays fell within this aperture, they 
were removed before flux extraction. In cases of crowded fields, aperture corrections were 
applied and ranged from 4" to 12" in diameter using the largest aperture possible without 
including contamination from neighboring sources. Uncertainties in the optical photometry 
were derived by estimating the internal night-to-night variations as well as the external errors 
(i.e. fits to the standard stars). A complete discussion of the error analysis can be found in 



^The central wavelengths for V, R, and / are 5475, 6425, and 8075A respectively. 



Henry et al.l (120041 ) . We adopt a total error of ±0.03 mag in each band. The final optical 
magnitudes are listed in Tabled] as well as the number of nights each object was observed. 



Infrared JHKs magnitudes and errors were extracted via Aladin from 2MASS and are 
also listed in Tabled! JHKs magnitude errors are, in most cases, significantly larger than for 
VRI, and the errors listed give a measure of the total photometric uncertainty (i.e. include 
both global and systematic components). In cases when the magnitude error is null, the star 
is near the magnitude limit of 2MASS and the photometry is not reliable. 



4. Analysis 
4.1. Modeling of Physical Parameters 



The pure hydrogen, pure helium, and mixed hy drogen and helium mo del atmospheres 
used to model the WDs are described at length in [Bergeron et al.l (120011 ) and references 
therein, while the heliuni -rich models appropriate for DQ and DZ stars are described in 
Dufour et al.l (120051 . l2007l ). respectively. The atmospheric parameters for each star are ob- 
tained by converting the optical VRI and infrared JHKs niagnitudes into observed fiuxes, 
and by comparing the resulting SEDs with those predicted from our model atmosphere cal- 
culations. The first step is accomplished by transfor ming the magii i tudes into average stellar 
fiuxes /j^ received at Earth using the calibration of iHolberg et al.l (120061 ) for photon count- 
ing devices. The observed and model fiuxes, which depend on T^s, logg, and atmospheric 
composition, are related by the equation 



/r = 47r {R/Df HI 



where R/D is the ratio of the radius of the star to its distance from Earth, and if™ is 
the Eddington fiux, properly averaged over the corresponding filter bandpass. Our fitting 



techn ique relies on the nonlinear least-squares method of Levenberg-Marquardt (jPress et al. 



19921 ). which is based on a steepest descent method. The value of is taken as the sum over 
all bandpasses of the difference between both sides of eq. (1), weighted by the corresponding 
photometric uncertainties. We consider only T^q and the solid angle to be free parameters, 
and the uncertainties of both parameters are obtained directly from the covariance matrix 
of the fit. In this study, we simply assume a value of log (7 = 8.0 for each star. 



As discussed in iBergeron et al.l (Il997l . I2OOII ). the main atmospheric constituent — hy- 
drogen or helium — is determined by comparing the fits obtained with both compositions, 
or by the presence of Ha in the optical spectra. For DQ and DZ stars, we rely on the 
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procedure outlined in iDufour et al.l (l2005l . 120071 ). respectively: we obtain a first estimate 
of the atmospheric parameters by fitting the energy distribution with an assumed value of 
the metal abundances. We then fit the optical spectrum to measure the metal abundances, 
and use these values to improve our atmospheric parameters from the energy distribution. 
This procedure is iterated until a self-consistent photometric and spectroscopic solution is 
achieved. 

The derived values for T^s for each object are listed in Table [TJ Also listed are the 
spectral types for each object determined based on their sp ectral features. The D As have 
been assigned a half- integer temperature index as defined by lMcCook &: Sioru (119991 ). where 
the temperature index equals 50,400/Tegf. As an external check, we compare in Figure [5] 
the photometric effective temperatures for the DA stars in Table 1 with those obtained by 
fitting the ob served Balmer l i ne prof iles (Figs. [2] and [3]) using t he spectroscopic tech nique 
developed by [Bergeron et al.l (Il992bl ) , and recently improved by iLiebert et al.l (120031 ) . Our 
grid of pure hydrogen, NLTE, and convective model atmospheres is also described in Liebert 
et al. The uncertainties of the spectroscopic technique are typically of 0.038 dex in log (7 
and 1.2% in T^s according to that study. We adopt a slightly larger uncertainty of 1.5% 
in Tefj (Spec) because of the problematic flux calibrations of the pre— 2006 data (see § 3.2). 
The agreement shown in Figure [5] is excellent, except perhaps at high temperatures where 
the photometric determinations become more uncertain. It is possible that the significantly 
elevated point in Figure El WD 0310—624 (labeled), is an unresolved double degenerate (see 
§ 4.2). We refrain here from using the log (7 determinations in our analysis because these 
are available only for the DA stars in our sample, and also because the spectra are not flux 
calibrated accurately enough for that purpose. 

Once the effective temperature and the atmospheric composition are determined, we 
calculate the a bsolu te visual magnitude of each star by combining the ne w calibration of 



Holberg et al.l (120061 ) with evolutionary models similar to those described in iFontaine et al. 
J200lh but with C/0 cores, g(He) = log Mrs/M^ = 10"^ ^nd g(H) = 10"^ (representative 
of hydrogen- atmosphere WDs), and q{B.e) = 10^^ and q{ii) = 10^^° (representative of 
helium-atmosphere WDs)!^!. By combining the absolute visual magnitude with the Johnson 
V magnitude, we derive a first estimate of the distance of each star (reported in Tabled]). 
Errors on the distance estimates incorporate the errors of the photometry values as well as 
an error of 0.25 dex in log g, which is the measured dispersion of the obs erved distribution 
using spectroscopic determinations (see Figure 9 of iBergeron et al.lll992bl ). 



Of the 33 new systems presented here, 5 have distance estimates within 25 pc. Four 



see http://www.astro.umontreal.ca/~bergeron/CoolingModels/ 
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more systems require additional attention because distance estimates are derived via other 
means. Three of these are hkely within 25 pc. All four are further discussed in the next 
section. In total, 20 WD systems (8 new and 12 known) are estimated (or determined) to be 
within 25 pc and one additional common proper motion binary system possibly lies within 
25 pc. 



4.2. Comments on Individual Systems 



Here we address unusual and interesting objects. 

WD 0121-429 is a DA WD that exhibits Zeeman splitting of Ha and H/3, thereby 
making its formal classification DAH. The SED fit to the photometry is superb, yielding 
a Teff of 6,369 ± 137 K. When we compare the strength of the absorption line trio with 
that predicted using the T^s from the SED fit, the depth of the absorption appears too 
shallow. Using the magnetic line fitting procedure outlined in iBergeron et al.l (jl992al ). we 
must include a 50% dilution factor to match the observed central line of Ha. In light of 
this, we utilized the trigonometric parallax distance determined via CTIOPI of 17.7 ± 0.7 
pc (Subasavage et al., in preparation) to further constrain this system. The resulting SED 
fit, with distance (hence luminosity) as a constraint rather than a variable, implies a mass of 
0.43 ± 0.03 Mp^ . Given the age of o ur Galaxy, the lowest mass WD that could have formed 
is ~0.47 Mq ( llben fc Renzinilll984j ). It is extremely unlikely that this WD formed through 
single star evolution. The most likely scenario is that this is a double degenerate binary 
with a magnetic DA component and a featureless DC comp onent (necessary to dilute the 
ab sorption at Ha), sirn ilar to G62-46 (IBergeron et al.lll993l ) and LHS 2273 (see Figure 33 



of iBergeron et al.l 119971 ). If this interpretation is correct, any number of component masses 
and luminosities can reproduce the SED fit. 

The spectrum and corresponding magnetic fit to the Ha lines (including the dilution) 
is shown in Figure O The viewing angle, i = 65°, is defined as the angle between the dipole 
axis and the line of sight (z = corresponds to a pole-on view). The best fit produces a 
dipole field strength, Bd = 9.5 MG, and a dipole offset, Qz = 0.06 (in units of stellar radius). 
The positive value of imphes that the offset is toward the observer. Only is moderately 
constrained, both i and can vary significantly yet still produce a reasonable fit to the data 
( IBergeron et al.lll992al ). 



WD 0310—624 is a DA WD that is one of the hottest in the new sample. Because of it's 
elevation significantly above the equal temperature line (solid) in Figure [5|, it is possible that 
it is an unresolved double degenerate with very different component effective temperatures. 
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In fact, this method has been used to identify unresolved double degenerate candidates (i.e. 
Bergeron et al.ll200ll ). 



WD 0511—415 is a DA WD (spectrum is plotted in Figure [2]) whose spectral fit 
produces a Tpff = 10,813 ±219 K and a log g = 8.21 ± 0.10 using the spectral fitting 



procedure of lLiebert et al.l (120031). This obj ect lies near the red edge of the ZZ Ceti instability 



strip as defined by iGianninas et al.l (120061 ). If variable, this object would help to constrain 
the cool edge of the instability strip in T^s, log g parameter space. Follow-up high speed 
photometry is necessary to confirm variability. 

WD 0622—329 is a DAB WD displaying the Balmer lines as well as weaker He i at 
4472 and 5876 A. The spectrum, shown in Figure (Tj is reproduced best with a model having 
Tgfj ~43,700 K. However, the predicted He ii absorption line at 4686 A for a WD of this 
Tgfj is not present in the spectrum. In contrast, the SED fit to the photometry implies a 
Teff of ~10,500 K (using either pure H or pure He models). Because the Tes values are 
vastly discrepant, we explored the possibility that this spectrum is not characterized by a 
single temperature. We modeled the spectrum assuming the object was an unresolved double 
degenerate. The best fit implies one component is a DB with Teff = 14,170 ± 1,228 K and 
the other component is a DA with Teff = 9,640 ± 303 K, similar to the unr esolved DA + 
DB degenerate binary PG 1115+166 analyzed by [Bergeron &: LiebertI (|2002| ). One can see 
from Figure [7] that the spectrum is well modeled under this assumption. We conclude this 
object is likely a distant (well beyond 25 pc) unresolved double degenerate. 

WD 0840-136 is a DZ WD whose spectrum shows both Ca ii (H & K) and Ca i (4226 
A) lines as shown in Figure [H Fits to the photometric data for different atmospheric com- 
positions indicate te mperatures of about 4800-5000 K. However, fits to the optical spectrum 
using the models of iDufour et al.l (120071 ) cannot repr oduce simultaneously all three calium 
lines. This problem is similar to that encountered by iDufour et al.l (120071 ) where the atmo- 
spheric parameters for the coolest DZ WDs were considered uncertain because of possible 
high atmospheric pressure effects. We utilize a photo metric relation rele vant for WDs of any 
atmospheric composition, which links My to {V — /) (ISalim et al.ll2004l ) to obtain a distance 
estimate of 19.3 ± 3.9 pc. 

WD 1054—226 was observed spectroscopically as part of the Edinbur gh-Cape (EC) 
blue object survey and assigned a spectral type of sdB-|- (IKilkenny et al.lll997l ). As is evident 
in Figure [HI the spectrum of this object is the noisiest of all the spectra presented here and 
perhaps a bit ambiguous. As an additional check, this object was recently observed using 
the ESO 3.6 m telescope and has been confirmed to be a cool DA WD (Bergeron, private 
communication) . 
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WD 1105—340 is a DA WD (spectrum is plotted in Figure [2]) with a common proper 
motion companion with separation of 30'.'6 at position angle 107.1°. The companion's spectral 
type is M4Ve with Vj = 15.04, Rkc = 13.68, Ikc = 11-96, J = 10.26, H = 9.70, and Ks 
= 9.41. In addition to the SED derived distance estimate for the WD, we utilize the main 



sequence distance relations of iHenry et al.l (120041 ) to estimate a distance to the red dwarf 



companion. We obtain a distance estimate of 19.1 ± 3.0 pc for the companion leaving 
open the possibility that this system may lie just within 25 pc. A trigonometric parallax 
determination is currently underway for confirmation. 

WD 1149—272 is the only DQ WD discovered in the new sample. This object was 
observed spectroscopically as part of the Edinburgh-Cape (EC) blue object s urvey for which 



no fe atures deeper than 5% were detected and was labeled a possible DC ([Kilkenny et al. 



19971 ). It is identified as having weak C2 swan band absorption at 4737 and 5165 A and is 
otherwise featureless. The DQ model reproduces the spectrum reliably and is overplotted in 
Figure m This object is characterized as having T^s = 6188 ± 194 K and a log (C/He) = 
-7.20 ± 0.16. 

WD 2008—600 is a DC WD (spectrum is plotted in Figure H]) that is flux deficient in 
the near infrared, as indicated by the 2MASS magnitudes. The SED fit to the photometry 
is a poor match to either the pure hydrogen or the pure helium models. A pure hydrogen 
model provides a slightly better match than a pure helium model, and yields a Tcs of ~3100 
K, thereby placing it in the relatively small sample of ultracool WDs. In order to discern the 
true nature of this object, we have constrained the model using the distance obtained from 
the CTIOPI trigonometric parallax of 17.1 ± 0.4 pc (Subasavage et al., in preparation). This 
object is then best modeled as having mostly helium with trace amounts of hydrogen (log 
(He/H) = 2.61) in its atmosphere and has a T^s = 5078 ± 221 K (see Figure [TOj) . A mixed 
hydrogen and helium composition is required to produce sufficient absorption in the infrared 
as a result of the collision-induced absorption by molecular hydrogen due to collisions with 
helium. Such mixed atmosph eric compositions hay e also been invoked to explain the infrared 



flux deficiency in LHS 112 6 (IBergeron et al.lll994h as well as SDSS 1337+00 and LHS 3250 



( iBergeron fc Leggettll2002h . While WD 2008-600 is likely not an ultracool WD, it is one of 
the brightest and nearest cool WDs known. Because the 2MASS magnitudes are not very 
reliable, we intend to obtain additional near-infrared photometry to better constrain the fit. 

WD 2138—332 is a DZ WD for which a calcium rich model reproduces the spectrum 
reliably. The spectrum and the overplotted fit are shown in the bottom panel of Figure 
El Clearly evident in the spectrum are the strong Ca 11 absorption at 3933 and 3968 A. A 
weaker Ca i line is seen at 4226A. Also seen are Mg i absorption lines at 3829, 3832, and 
3838 A (blended) as well as Mg i at 5167, 5173, and 5184 A (also blended). Several weak Fe i 
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lines from 4000A to 4500A and again from 5200A to 5500A are also present. The divergence 
of the spectrum from the fit toward the red end is likely due to an imperfect flux calibration 
of the spectrum. This object is characterized as having Teg = 7188 ± 291 K and a log 
(C a/He) = —8.64 ± O-IS. T he metallicity ratios are, at first, assumed to be solar (as defined 



by iGrevesse fc Sauvallll998l ) and, in this case, the quality of the fit was sufficient without 
deviation. The corresponding log (Mg/He) = —7.42 ± 0.16 and log (Fe/He) = —7.50 ± 0.16 
for this object. 

WD 2157—574 is a DA WD (spectrum is plotted in Figure[3]) unique to the new sample 
in that it displays weak Ca ii absorption at 3933 and 3968 A (H and K) thereby making its 
formal classification a DAZ. Possible scenarios that enrich the atmosph eres of DAZs inc lude 



accretion via (1) debris disks, (2) ISM, and (3) cometary impacts (see lKilic et al.ll2006l and 
references therein). The 2MASS Ks magnitude is near the faint limit and is unreliable, but 
even considering the J and H magnitudes, there appears to be no appreciable near-infrared 
excess. While this may tentatively rule out the possibility of a debris disk, this object would 
be an excellent candidate for far-infrared spaced-based studies to ascertain the origin of the 
enrichment. 



5. Discussion 

WDs represent the end state for stars less massive than ~8 and are therefore rel- 
atively numerous. Because of their intrinsic faintness, only the nearby WD population can 
be easily characterized and provides the bench mark upon which WD stellar astrophysics is 



based. It is clear from this work and others (e.g. iHolberg et al.ll2002l : iKawka fc Vennesll2006l ) 
that the WD sample is complete, at best, to only 13 pc. Spectroscopic confirmation of new 
WDs as well as trigonometric parallax determinations for both new and known WDs will 
lead to a more complete sample and will push the boundary of completeness outward. We 
estimate that 8 new WDs and an additional 12 known WDs without trigonometric parallaxes 
are nearer than 25 pc, including one within 10 pc (WD 0141—675). Parallax measurements 
via CTIOPI are underway for these 20 objects to confirm proximity. This total of 20 WDs 
within 25 pc constitutes an 18% increase to the 109 WDs with trigonometric parallaxes > 
40 mas. 

Evaluating the proper motions of the new and known samples within 25 pc indicates 
that almost double the number of systems have been found with /x < I'.'O yr^^ than with 
/i > I'.'O yr^^ (13 vs 7, see Table [2]). The only WD estimated to be within 10 pc has /i > 
I'.'O yr-\ ahhough WD 1202-232 is estimated to be 10.2 ± 1.7 pc and it's proper motion is 
small (/i = 0'.'227 yr'^). 
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Because this effort focuses mainly on the southern hemisphere, it is hkely that there 
is a significant fraction of nearby WDs in the northern hemisphe re that have also gon e 



undetected. With the recent release of the LSPM-North Catalog (iLepine fc Shard l2005l ) . 
these objects are identifiable by employing the same techniques used in this work. The 
challenge is the need for a large scale parallax survey focusing on WDs to confirm proximity. 
Since the HIPPARCOS mission, only s ix WD trigonometric parallaxes have been published 



( iHambly et al.lll999l : ISmart et al.ll2003[ ). and of those, only two are within 25 pc. The USNO 
parallax program is in the process of publishing trigonometric parallaxes for ~130 WDs, 
mostly in the northern hemisphere, although proximity was not a primary motivation for 
target selection (Dahn, private communication). 

In addition to further completing the nearby WD census, the wealth of observational 
data available from this effort provides reliable constraints on their physical parameters 
(i.e. Teff, log mass, and radius). Unusual objects are then revealed, such as those dis- 
cussed in § 4.2. In particular, trigonometric parallaxes help identify WDs that are overlu- 
minous, as is the case for WD 0121—429. This object, and others similar to it, are excellent 
candidates to provide insight into binary evolution. If they can be resolved using high res- 
olution astrometric techniques (i.e. speckle, adaptive optics, or interferometry via Hubble 
Space Telescope's Fine Guidance Sensors), they may provide astrometric masses, which are 
fundamental calibrators for stellar structure theory and for the reliability of the theoretical 
WD mass-radius and initial-to-final-mas s relationships. To da te, only four WD astrometric 



masses are known to better than ~ 5% (jProvencal et al.lll998l ) 



One avenue that is completely unexplored to date is a careful high resolution search 
for planets around WDs. Theory dictates that the Sun will become a WD, and when it 
does, the outer planets will remain in orbit (not without transformations of their own, of 
course). In this scenario, the Sun will have lost more than half of its mass, thereby amplifying 
the signature induced by the planets. Presumably, this has already occurred in the Milky 
Way and systems such as these merely await detection. Because of the faintness and spectral 
signatures of WDs (i.e. few, if any, broad absorption lines), current radial velocity techniques 
are inadequate for planet detection, leaving astrometric techniques as the only viable option. 
For a given system, the astrometric signature is inversely related to distance (i.e. the nearer 
the system, the larger the astrometric signature). This effort aims to provide a complete 
census of nearby WDs that can be probed for these astrometric signatures using future 
astrometric efforts. 
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A. Appendix 

In order to ensure correct cross-referencing of names for the new and known WD systems 
presented here, Table[3]lists additional names found in the literature. Objects for which there 
is an NLTT designation will also have the corresponding L or LP designations found in the 
NLTT catalog. This is necessary because the NLTT designations were not published in the 
original catalog, but rather are the record numbers in the electronic version of the catalog 
and have been adopted out of necessity. 
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Table 1. Optical and Infrared Photometry, and Derived Parameters for New and Known White Dwarfs. 
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Table 1 — Continued 



WD Vj Rc Ic # CTj H CTH Ks CTKg Tett Comp Dist SpT 

Name Obs ' (K) (pc) 



1814+134 15.85 15.34 14.86 2 14.38 0.04 14.10 0.06 14.07 0.06 5313± 115 H 15.6±2.5 DA9.5 

2040-392 13.74 13.77 13.68 2 13.77 0.02 13.82 0.03 13.81 0.05 10811± 325 H 23.1±4.0 DA4.5 

2211-392 15.91 15.61 15.24 2 14.89 0.03 14.64 0.05 14.56 0.08 6243± 167 H 23.5±4.0 DAS.O 

2226-754A 16.57 15.93 15.33 2 14.66 0.04 14.66 0.06 14.44 0.08 4230± 104 H 12.8±2.0 DC 

2226-754B 16.88 16.17 15.51 2 14.86 0.04 14.82 0.06 14.72 0.12 4177± 112 H 14.0±2.2 DC 



Known White Dwarfs without a Trigonometric Parallax Estimated to be Beyond 25 pc 
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^Distance via SED fit (not listed) is underestimated beeause object is likely an unresolved double degenerate with one magnetic component (see § 4.2). Instead, we 
adopt the trigonometric parallax distance of 17.7 i 0.7 pc derived via CTIOPI. 

'^Distance via SED fit (not listed) is underestimated because object is likely a distant (well beyond 25 pc) unresolved double degenerate (see § 4.2). 

'^Distance via SED fit (not fisted) is underestimated because object is likely a distant (weU beyond 25 pc) unresolved double degenerate with components of type DA 
and DB (see § 4.2). Temperatures derived from the spectroscopic fit yield 9,640 it 303 K and 14,170 d: 1,228 K for the DA and DB respectively. 

'^Object is likely cooler than T^^ff ^^^ 5000 K and t he th eoretical models do not provide an accurate treatment at these temperatures (sec § 4.2). Instead, we use the 
linear photometric distance relation of lSalim et al"] j2004| ) and obtain a distance estimate of 19.3 zb 3.9 pc. 

'^This object was observed as part of the Edinburgh-Cape survey and was classified as a sdB-|- jKilkennv et al.lfl997l l. 

^Distance of 19.1 zb 3.0 pc is estimated using VRIJHKs for the common proper motion companion M dwarf and the relations of iHenrv et aT ] i2004l 1. Syst em IS 
possibly within 25 pc. (see § 4.2). 

^Distance estimate is undetermined. Instead, we adopt the distance measured via trogonometric parallax of 17.1 it 0.4 pc (see § 4.2). 

"^Not listed in lMcCook 8z SioTil i 19991 ) but identified as a DC/DQ WD bv lHenrv et al.| i2002fl . We obtained blue spectra that show Ca II H &; K absorption and classify 
this object as a DZ. 

^The SED fit to the photometry is marginal. This object displays deep swan band absorption that significantly affects its measured magnitudes. 
J Not listed in lMcCook &€ SionI il999D but identified as a WD by |Luyten| | |1949| ). Spectral type is derived from our spectra. 

^As of mid-2004, object has moved onto a background source. Photometry is probably contaminated, which is consistent with the poor SED fit for this object. 
^Spectral type was determined using spectra published bv lSchoIz et aT] j2002| ) . 
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Table 2. Distance Estimate Statistics for New and Known White Dwarfs. 



Proper motion 


d < 10 pc 


10 pc < d < 25 pc 


d > 25 pc 


II > I'.'O yr-i 


1 


6 


1 


I'/O yr-i > > <y.'8 yr-i 











Of.'S yr-i > fi> O'.'e yr-i 





2 


2 


O'.'e yr-i > fi> 0'.'4 yr"! 





6 


11 


0'.'4 yr-i > ii> O'.'IS yr"!.... 





5 


22 


Total 


1 


19 


36 
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Table 3. Astrometry and Alternate Designations for New and Known White Dwarfs. 



RA 
(J2000.0) 



Dec 
(J2000.0) 



PM 
{arcsec yr~ 



PA 
(deg) 



Ref 



Alternate Names 



New Spectroscopically Confirmed White Dwarfs 



0034- 


602 


00 


36 


22 


31 


-59 


55 


27 


5 





280 


069 





L 


NLTT 1993 


_ 


LP 122-4 


_ 




0121- 


429 


01 


24 


03 


98 


-42 


40 


38 


5 





538 


155 


2 


L 


LHS 1243 


_ 


NLTT 4684 


_ 


LP 991-16 


0216- 


398 


02 


18 


31 


51 


—39 


36 


33 


2 





500 


078 


6 


L 


LHS 1385 


_ 


NLTT 7640 


_ 


LP 992-99 


0253- 


755 


02 


52 


45 


64 


-75 


22 


44 


5 





496 


063 


5 


S 


SCR 0252-7522 


_ 




_ 




0310- 


624 


03 


11 


21 


34 


-62 


15 


15 


7 





416 


083 


3 


S 


SCR 0311-6215 


= 




= 




0344+ 


014 


03 


47 


06 


82 


+ 01 


38 


47 


5 





473 


150 


4 


s 


LHS 5084 




NLTT 11839 




LP 593-56 


0404- 


510 


04 


05 


32 


86 


-50 


55 


57 


8 





320 


090 


7 


p 


LEHPM 1-3634 










0601- 


665 


05 


02 


43 


43 


-65 


26 


35 


2 





280 


191 


9 


p 


LBHPM 1-3866 










0611- 


415 


05 


13 


27 


80 


-41 


27 


51 


7 





292 


004 


4 


p 


LEHPM 2-1180 










0525- 


311 


05 


27 


24 


33 


-31 


06 


55 


7 





379 


200 


7 


p 


NLTT 15117 




LP 892-45 




LEHPM 2-521 


0607- 


530 


06 


08 


43 


81 


-53 


01 


34 


1 





246 


327 


6 


p 


LEHPM 2-2008 










0622- 


329 


06 


24 


25 


78 


-32 


57 


27 


4 





187 


177 


7 


p 


LEHPM 2-5035 










0821- 


669 


08 


21 


26 


70 


-67 


03 


20 


1 





758 


327 


6 


s 


SCR 0821-6703 










0840- 


136 


08 


42 


48 


45 


-13 


47 


13 


1 





272 


263 





s 


NLTT 20107 




LP 726-1 






1016- 


308 


10 


18 


39 


84 


-31 


08 


02 








212 


304 





L 


NLTT 23992 




LP 904-3 




LEHPM 2-6779 


1054- 


226 


10 


56 


38 


64 


-22 


52 


55 


9 





277 


349 


7 


P 


NLTT 25792 




LP 849-31 




LEHPM 2-1372 


1105- 


340 


11 


07 


47 


89 


-34 


20 


51 


4 





287 


168 





s 


SCR 1107-3420A 










1149- 


272 


11 


51 


36 


10 


-27 


32 


21 








199 


278 


3 


p 


LEHPM 2-4051 










1243- 


123 


12 


46 


00 


69 


-12 


36 


19 


9 





406 


305 


4 


s 


SCR 1246-1236 










1316- 


215 


13 


19 


24 


72 


-21 


47 


55 








467 


179 


2 


s 


NLTT 33669 




LP 854-50 




WT 2034 


1436- 


781 


14 


42 


51 


54 


-78 


23 


53 


6 





409 


272 





s 


NLTT 38003 




LP 40-109 




LTT 5814 


1462- 


310 


14 


55 


23 


47 


-31 


17 


06 


4 





199 


174 


2 


p 


LEHPM 2-4029 










1647- 


327 


16 


50 


44 


32 


-32 


49 


23 


2 





526 


193 


8 


L 


LHS 3245 




NLTT 43628 




LP 919-1 


1742- 


722 


17 


48 


31 


21 


-72 


17 


18 


5 





294 


228 


2 


P 


LEHPM 2-1166 










1946- 


273 


19 


49 


19 


78 


-27 


12 


25 


7 





213 


162 





L 


NLTT 48270 




LP 925-53 






2008- 


600 


20 


12 


31 


75 


-59 


56 


51 




1 


440 


165 


6 


S 


SCR 2012-5956 










2008- 


799 


20 


16 


49 


66 


-79 


45 


53 








434 


128 


4 


S 


SCR 2016-7945 










2035- 


369 


20 


38 


41 


42 


-36 


49 


13 


5 





230 


104 





L 


NLTT 49589 




L 495-42 




LEHPM 2-3290 


2103- 


397 


21 


06 


32 


01 


-39 


35 


56 


7 





266 


151 


7 


P 


LEHPM 2-1671 










2138- 


332 


21 


41 


57 


56 


-33 


00 


29 


8 





210 


228 


5 


P 


NLTT 51844 




L 570-26 




LEHPM 2-3327 


2157- 


574 


22 


00 


45 


37 


-57 


11 


23 


4 





233 


252 





P 


LEHPM 1-4327 










2218- 


416 


22 


21 


25 


37 


-41 


25 


27 








210 


143 


4 


P 


LEHPM 1-4598 










2231- 


387 


22 


33 


54 


47 


-38 


32 


36 


9 





370 


220 


5 


P 


NLTT 54169 




LP 1033-28 




LEHPM 1-4859 



Known White Dwarfs without a Trigonometric Parallax Estimated to be Within 25 pc 



0141-676 ... 
0806-661 ... 
1009-184 ... 
1036-204 ... 
1202-232 ... 
1315-781 ... 
1339-340 ... 
1766+143 ... 
1814+134 ... 
2040-392 ... 
2211-392 ... 
2226 -764 A.. 
2226-764B.. 



0024-556.. 
0150+256.. 
0255-705.. 
0442-304.. 
0928-713.. 
1143-013.. 



01 


43 


00 


98 


-67 


18 


30 


3 


1 


048 


197.8 


L 


LHS 145 




NLTT 5777 




L 88-59 


08 


06 


53 


76 


-66 


18 


16 


6 





464 


131.4 


S 


NLTT 19008 




L 97-3 






10 


12 


01 


88 


-18 


43 


33 


2 





519 


268.2 


S 


WT 1759 




LEHPM 2-220 






10 


38 


55 


57 


-20 


40 


56 


7 





628 


330.3 


L 


LHS 2293 




NLTT 24944 




LP 790-29 


12 


05 


26 


66 


-23 


33 


12 


1 





227 


002.0 


L 


NLTT 29555 




LP 852-7 




LEHPM 2-1894 


13 


19 


25 


63 


-78 


23 


28 


3 





477 


139.2 


S 


NLTT 33551 




L 40-116 






13 


42 


02 


88 


-34 


15 


19 


4 


2 


547 


296.7 


Le 


PM J13420-3415 










17 


58 


22 


90 


+ 14 


17 


37 


8 


1 


014 


235.4 


Le 


LSR 1768+1417 










18 


17 


06 


48 


+13 


28 


25 





1 


207 


201.6 


Le 


LSR 1817+1328 










20 


43 


49 


21 


-39 


03 


18 








306 


179.0 


L 


NLTT 49752 




L 495-82 






22 


14 


34 


75 


-38 


59 


07 


3 


1 


066 


110.1 


O 


WD J2214-390 




LEHPM 1-4466 






22 


30 


40 


00 


-76 


13 


55 


3 


1 


868 


167.5 


s 


SSSPM J2231-7514 










22 


30 


33 


55 


-76 


15 


24 


2 


1 


868 


167.5 


s 


SSSPM J2231-7515 














Known White Dwarfs without a Trigonometric 


Parallax 


Estimated to be Beyond 25 pc 






00 


26 40.69 


-55 


24 44.1 


0.580 


211.8 


L 


LHS 1076 




NLTT 1415 




L 170-27 


01 


52 


51.93 


+25 53 40.7 


0.220 


076.0 


L 


NLTT 6275 




G 94-21 






02 


56 


17.22 


-70 22 


10 


8 


0.682 


097.9 


L 


LHS 1474 




NLTT 9485 




L 54-5 


04 44 29.38 


-30 21 


14.2 


0.196 


199.5 


P 


NLTT 13882 




LP 891-65 




HE 0442-3027 


09 


29 


07.97 


-71 


33 


58 


8 


0.439 


320.2 


s 


NLTT 21957 




L 64-40 






11 


46 


25.77 


-01 


36 36 


8 


0.563 


140.2 


s 


LHS 2455 




NLTT 28493 
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Table 3 — Continued 



WD Name RA Dec PM PA Rcf Alternate Names 

(J2000.0) (J2000.0) (arcsec yr"!) (deg) 

1237-230 12 40 24.18 -23 17 43.8 1.102 219.9 L LHS 339 = NLTT 31473 = LP 853-15 

1314-153 13 16 43.59 -15 35 58.3 0.708 196.7 L LHS 2712 = NLTT 33503 = LP 737-47 

1418-088 14 20 54.93 -09 05 08.7 0.480 266.8 S LHS 5270 = NLTT 37026 = ■ ■ ■ 

1447-190 14 50 11.93 -19 14 08.7 0.253 285.4 P NLTT 38499 = LP 801-14 = LEHPM 2-1835 

1607-250 16 10 50.21 -25 13 16.0 0.209 314.0 L NLTT 42153 = LP 861-31 = ■ ■ ■ 



Refere nces. — (L) lLuvtenlll979ai rbl (Lei iLepine et aLll2003l . iLepine et aLll2005l . (O) lOppenheimer et al1l200ll . (P) iPokcrnv et aLl 

120041 . (S) ISubasavage et aLll2005al lhl this work 
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Fig. 1. — Reduced proper motion diagram used to select WD candidates for spectrosco pic 



follow-up. Plotted are the new high proper motion objects from lSubasavage et al.l (l2005a| ]bl). 
The line is a somewhat arbitrary boundary between the WDs (below) and the subdwarfs 
(just above). Main sequence dwarfs fall above and to the right of the subdwarfs, although 
there is significant overlap. Asterisks indicate the 33 new WDs reported here. Three dots 
in the WD region are deferred to a future paper. The point labeled "sd" is a confirmed 
subdwarf contaminant of the WD sample. 

Fig. 2. — Spectral plots of the hot (Tefj > 10000 K) DA WDs from the new sample, plotted 
in descending T^^ as derived from the SED fits to the photometry. Note that some of the 
flux calibrations are not perfect, in particular, at the blue end. 

Fig. 3. — Spectral plots of cool (Tcfj < 10000 K) DA WDs from the new sample, plotted in 
descending Teg as derived from the SED fits to the photometry. Note that some of the flux 
calibrations are not perfect, in particular, at the blue end. 

Fig. 4. — Spectral plots of the four featureless DC white dwarfs from the new sample, plotted 
in descending Tgfj as derived from the SED fits to the photometry. Note that some of the 
flux calibrations are not perfect, in particular, at the blue end. 

Fig. 5. — Comparison plot of the values of Tgfj derived from photometric SED fitting vs 
those derived from spectral fitting for 25 of the DA WDs in the new sample. The solid line 
represents equal temperatures. The elevated point, 0310—624, is discussed in § 4.2. 

Fig. 6. — Spectral plot of WD 0121—429. The inset plot displays t he spectrum (light line ) 



in the Ha region to which a magnetic fit (heavy line), as outlined in lBergeron et al.l (jl992af ). 
was performed using the Tgfj obtained from the SED fit to the photometry. The resulting 
magnetic parameters are listed below the fit. 

Fig. 7. — Spectral plot of WD 0622—329. The inset plot displays the spectrum (light line) 
in the region to which the model (heavy line) was fit assuming the spectrum is a convolution 
of a DB component and a slightly cooler DA component. Best fit physical parameters are 
listed below the fit for each component. 

Fig. 8. — {top panel) Spectral plot of WD 0840—136. The DZ model failed to reproduce 
the spectrum presumably because this object is cooler than T^q ~ 5000 K where additional 
pressure effects, not included in the model, become important, {bottom panel) Spectral plot 
of WD 2138—332. The inset plot displays the spectrum (light line) in the region to which 
the model (heavy line) was fit. 
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Fig. 9. — Spectral plot of WD 1149—272. The inset plot displays the spectrum (light hne) 
in the region to which the model (heavy line) was fit. 

Fig. 10. — Spectral energy distribution plot of WD 2008—600 with the distance constrained 
by the trigonometric distance of 17.1 ± 0.4 pc. Best fit physical parameters are hsted below 
the fit. Points are fit values; error bars are derived from the uncertainties in the magnitudes 
and the parallax. 
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; log (C/He) = -7.20 ± 0.16 ! 
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